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Abstract 

We consider the process e + e~ — ► N1N2, where N% and N 2 are heavy Ma- 
jorana particles, with relative CP given by rjcp = +1 or —1, decaying sub- 
sequently via Ni,N 2 -+ W ± e T . We derive the energy and angle correlation 
of the dilepton final state, both for like-sign (e^e^) and unlike-sign (e~e + ) 
configurations. Interesting differences are found between the cases rjcp = +1 
and —1. The characteristics of unlike-sign e + e~ dileptons originating from 
a Majorana pair N1N2 are contrasted with those arising from the reaction 
e + e~ — > -/ViV — > W + e~W~e + , where A^A^ is a Dirac particle-antiparticle 
pair. 
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I. INTRODUCTION 



In an interesting paper [||], Kogo and Tsai have analysed the reaction e + e~ — *> N1N2, 
where Ni j2 are heavy Majorana neutrinos, and compared the cases where the relative CP 
of Ni and N 2 is rjcp = +1 and —1. It was found that the two cases differ in threshold 
behaviour, in angular distribution, and in the dependence on the spin-directions of Ni and 
N 2 . A comparison was also made between the Majorana process and the Dirac process 
e + e~ — > NN, where NN is a Dirac particle-antiparticle pair. A related analysis was carried 
out in Ref. [fj. (The contrast between Majorana neutrinos and Dirac neutrinos has been 
the subject of several other papers (e.g. 0-§J) and monographs ( 0H)) 

In the present paper, we examine how the differences between the cases rjcp = +1 and 
— 1 propagate to the decay products of N\ and N 2 , assuming the decays to take place via 
Ni t 2 — > W ± e T . We focus on the like-sign lepton pair created in the reaction chain e + e~ — > 
N\N 2 — > W + e~W + e~ , which is a characteristic signature of Majorana pair production. We 
derive, in particular, the correlation in the energies of the e~e~ pair, and in their angles 
relative to the e + e~ axis. Interesting differences are found between the cases rjcp = +1 
and — 1. We also examine the behaviour of the unlike-sign dileptons e + e~, comparing the 
Majorana cases with dileptons created in the production and decay of a Dirac NN pair, i.e. 
e+e" -> NN -> W + e~W~e + . 



II. CHARACTERISTICS OF THE REACTION e+e 



N t N 2 



The analysis of Ref. [I] was carried out in the context of the simple production mechanism 
for e + e~ — > N±N 2 shown in Fig. 1, and we begin by recapitulating the essential results. The 
interaction Lagrangian is taken to be 
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+ a N N 2 (x) 7m _(l - 75 ) JVi(x) 



ZHx) 



(2.1) 



where cy, ca and a at may be regarded as real phenomenological parameters. (For the 
standard Z-boson, cy = — 1/2 + 2 sin 2 6^, ca = —1/2). The matrix element for Majorana 
neutrinos (with momenta and spins as indicated in Fig. 1) is 



M m = -ia N {^^-? 3t Af 
2 cos u w 



UtAqi) Iv ^(l-7s) v t2 {q 2 ) A 2 
- u t2 {q 2 ) lu ^(1 - 75) v tl (qi) Ai 



(2.2) 



where 



3l = v S2 (p 2 ) In (c v - ca1 5 ) u sl (pi) (2.3) 



and 



z q 2 - m 2 z +im z T z " 



Assuming CP-invariance the factors Ai, A2 in Eq. (2.2) are such that A1A2 = +1(— 1) when 
Ni and N 2 have the same (opposite) CP-parity ||. Rewriting the second term in Eq. (2.2) 
as 

WfaCftj) lv -(1 -75) = u tl (qi) j„ -(1 + j 5 ) v t2 (q 2 ) , (2.5) 

we observe that the current of the Majorana neutrinos is pure axial vector when N% and N 2 
have the same CP-parity (r]cp = \±\ 2 = +1), and pure vector when they have opposite CP 
(Vcp — A1A2 = —1). In comparison, the matrix element for the Dirac process e + e~ — > NN 
is 

M d = -ia N { T J ~ ir Y f» Af u tl {qi) lu ~(1 - 75) v t2 (q 2 ) . (2.6) 
2 cos t>w 2 

The differential cross section for e + e~ — > NiN 2 , for general masses mi and m 2 , and for 
arbitrary polarizations n and n' of the two neutrinos is given in the Appendix. In Sec. 5 we 
compare our formulas with those of Ref . [|1] , and with special cases treated in other papers. 
Here we specialise to the case m\ = m 2 = m^v, for which the cross-section (da/dQ) in the 
cases rjcp = +1 and —1 is 



{-£)+ = \o* P 3 { fi [{n y n' y - n x n' x )S 2 + (1 + n z n' z )(l + C 2 )] - f 2 2{n z + n' z )C } , (2.7) 

=cr p[f 1 [2-p 2 + C 2 f3 2 + n z n' z ([3 2 + C 2 (l/ 7 2 + 1)) + n,<5 2 (l/ 7 2 + 1) 

-n y n' y S 2 (3 2 + (n x n' z + n' x n z )2SC h 2 } - f 2 [2(n x + n' x )Sh 2 + 2(n z + n' z )C] } . (2.8) 

For comparison, the differential cross section of the Dirac process e + e~ — * NN is 

{%)* = \°*P{ h [(1 + C 2 i3 2 ) - (n z + n'Ml + C 2 ) (2.9) 
-{n x + ri x )SCPh + n z n' z {C 2 + (3 2 ) + {n x n' z + n x n' x )SC h + n x n' x S 2 / 7 2 ] + 
/ 2 [2C/3 - (n 2 + n'JC(l + /? 2 ) - (n x + <)S/ 7 + 2n z n' z Cf3 + (n x < + n' x n z )SPh\ }. 

The symbols in Eqs. (2.7)-(2.9) are defined as follows: 

s , ) 9 = (l-4my a ) 1 / 2 , 7 = (l-/3 2 )- 1/2 

(2.10) 

C = cos6> , S^sinfl , f 1 = 2(c 2 v + c 2 A ) , f 2 = 4 c v c A , 

# being the scattering angle of Ni (or JV) with respect to the initial e~ direction. The 
co-ordinate axes are defined so that the momentum- and spin-vectors of Ni and N 2 in the 
e + e~ cm. frame have the components 

= ( 7 m, 0, 0, 7/3m) , # = {~f(3n z ,n x ,n y ,~in z ) , (2.11) 
?2 = ( 7 m, 0, 0, - 7 /3m) , f£ = (- 7 /3<> <, nj,, 7 <) . 

Inspection of Eqs. (2.7)-(2.9) reveals several interesting features: 

(a) The Majorana cases '+' and '— ' have different dependence on the spin-vectors n and 
ft', and different angular distributions, even after the spins n and n' are summed 
over. These differences stem from the fact that the matrix element M. m in Eq. (2.2) 
effectively involves an axial vector current iV l7jti75 7V2 when A1A2 = +1 and a vector 
current N\^ lJl N 2 when A1A2 = — 1. 
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cr = 



r 2 n, 2 

5127T 2 



s — m z + im z T z 



(b) The Majorana cases '+' and '— ' differ from the Dirac case 'd', in which the current of 
the AA-pair has a V-A structure A7 M |(1 — 75) A. This difference persists even if the 
spins of the heavy neutrinos are summed over, in which case 



E(^) + = 2a /3 3 [ /l (l + C 2 ) 



A a, (3 / 1 (2-/3 2 + C 2 /3 2 ) 



E(— ) 

n,n' 



(2.12) 



Whereas the spin-averaged Majorana cross sections are forward-backward symmetric, 
the Dirac process has a term linear in cos 9, with a coefficient proportional to j 2 = 
4 cyCA- Eq. (2.12) also shows that the threshold behaviour is (3 3 , (3 and (3 for the cases 
'+', '— ' and 'd' respectively. In the asymptotic limit (3 — > 1 the Majorana cases '+' 
and '— ' have the same angular distribution (1 + C 2 ), distinct from that of the Dirac 
process. 

(c) In the high energy limit f3 — > 1, the Dirac process e + e~ — > iViV has a spin-dependence 
given by 



;^) d = ia /3[l-K + <)+nX][/i(l + C 2 ) + 2f 2 C 



(2.13) 



The fact that only the longitudinal components (n z and n' z ) of the N, N spins appear in 
this expression is consistent with the expectation that relativistic Dirac neutrinos are 
eigenstates of helicity. The fact that the cross section (2.13) vanishes when n z = —1, 
n' z — +1 confirms the expectation that for a V-A current the and N are produced 
in left-handed and right-handed states, respectively. By comparison, the Majorana 
processes e + e~ — > N ± N 2 , for rjcp = +1 and —1, have the high energy behaviour 

(0-1) 



, c?cr 



)+ = { /i [(1 + C 2 )(l + n 2 <) + S 2 (n y n' y - n x n' x )} - 2f 2 C (n z + n' z ) } , (2.14) 
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(^)- = { h [(1 + C 2 )(l + n z n' z ) + S\n x n' x - n y n' y )} - 2f 2 C (n z + n' z ) } . (2.15) 



Contrary to the Dirac case, the Majorana reactions have an explicit dependence on n x , 
n y and n' x , n' y , reflecting the fact that a relativistic Majorana particle with m N ^ is 
not necessarily an eigenstate of helicity, and can have a spin pointing in an arbitrary 
direction. The Majorana cases '+' and '— ' differ in the sign of the term proportional 
to S 2 , which contains the transverse (x- and y-) components of the neutrino spins. It 
is with the purpose of exposing the subtle differences in the spin state of the NiN 2 and 
NN systems that we investigate in the following sections the dilepton final state created 
by the decays of the heavy neutrinos via N lj2 — > W ± e T and N(N) — > W + e~(W~e + ). 

III. LIKE-SIGN DILEPTONS: THE REACTION e+e" -► iViiVa -► W + W + e~e~ 

As seen in the preceding, the spin state and the angular distribution of the Majorana 
pair produced in e + e~ — > NiN 2 depends on the relative CP-parity rjcp of the two particles. 
We wish to see how these differences manifest themselves in the decay products of Ni and 
N 2 . To this end, we assume that m N > m w , and that the simplest decay mechanism is 
N lj2 — > W z ^e ± . In particular, the reaction sequence e + e~ NxN 2 — > W + W + e~ e~ leads 
to the appearance of two like-sign leptons in the final state, an unmistakable signature of 
Majorana pair production. (For the purpose of this paper we assume that the IT^-bosons 
decay into quark jets, thus avoiding the complications of final states with 3 or 4 charged 
leptons.) 

We have calculated the amplitude of the process e + e~ — > N\N 2 — > W + W + e~e~ , depicted 
in Fig. 2, assuming a decay interaction (a' N and a" N beeing real parameters) 



C 2 {x) 
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a N e{x) llx -{l- lb )N 1 {x) W^(x) 

+ d N N 1 {x) 1 » l -{l- lb )e{x) W^(x) 
+ a N e{x) 1 ^{l- lb )N 2 {x) W^{x) 



V2 
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+ a" N N 2 (x) lfl ^(l- l5 )e(x) W» + (x) 



This amplitude has the form (see Appendix for details) 



(3.1) 



M=iA j; A£ 



q( — m\ + imiYi q 2 ~~ m 2 + ^ m 2^2 
m 2 \ 2 u tl (h) 7 P fa 7^ 7 CT ^(1 + 75) v t2 (k 2 ) 

m 1 \ 1 u t2 (k 2 ) 7cr ^2 lv 1 P tt(1 + 75) ^(fci) 



2 



^3)^4) , (3.2) 



where 



A = a " a " a *'w^- (3 - 3) 

Using the narrow-width approximation for the Ni, N 2 propagators, and specializing to the 
case mi = m 2 = to at, we obtain the following expression for the squared matrix element 
(summed over final and averaged over initial spins), the subscript in M.± denoting r\ C p = ±1 
(q = Pi +P2, I = Pi -P2): 

\ A \ 2 1 T xf 2 2^ f r/2 2 x m2 



/1 • ( T(m 2 N - m^) 2 (m^ + 2m 2 w f ■ s T 4(m^ - 2m 2 w ) 



4 

//t 

2 \2 



s {kik 2 ){qiq2) ~ s (k 1 q 2 )(k2qi) - (k 1 k 2 )(qiq)(q2q) + (hh) (qj) (q 2 l) 
+ (hq2)(hq)(qiq) - (k 1 q 2 )(k 2 l)(qil) + {kiq){k 2 q 1 ){q 2 q) - (hi) (k 2 q^ (q 2 l) 
~ (hq)(k 2 q)(qiq 2 ) + (hi) (k 2 l) ( qi q 2 ) ± m 2 N ((hq)(k 2 q) - (hl)(k 2 l)) 
+ 2 (m 2 N - m 2 w )(m 2 N - 2m 2 v ) 2 (hq)(q 2 q) - (hl)(q 2 l) + (hq)(qiq) - (k 2 l)(qil) 

(92?) - (9i0(92/) 
-2/2 • (m 2 N - m 2 w )(m A N - 4to^) ^ ±(hq)(qj) - (hq)(q 2 l) 
T(hl)(qiq) + (hl)(q 2 q) ± (hq)(q 2 l) ~ (hq)(qil) 

T(M)(929) + (M)(9i9) ) } • (3.4) 
If the final state is e + e + instead of e~e~, we replace f'2 — > — /b in the above equation. 



1 o 2/2 2 \2 

+ 8 m w (m N - m w ) 
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The expression for |.M±| 2 can be integrated over the phase space of W + and W + (i.e. 
over the momenta ks(= qi — k\) and k 4 (= q 2 — k 2 ), in order to obtain the spectra in the 
lepton variables ki and k 2 . Defining the four-vectors k\ and k 2 in the e + e~ cm. frame by 

ki = Ei (1, sin 6\ cos 0i, sin Q\ sin 0i, cos d\) , 

k 2 = E 2 (1, sin 9 2 cos 4> 2) sin 9 2 sin 2 , cos 9 2 ) , (3.5) 

we have been able to derive the correlated distribution of the energies E\ and E 2) as well as 
the correlation of the variables cos Q\ and cos 9 2 measured relative to the e~ beam direction. 



A. Energy Correlation 

The normalized spectrum in the energies of the dilepton pair e~e~ is (£\ >2 = Ei^/m^) 

-•(-77^) = A/" [a + b(£i + £ 2 ) + c(£i + £ 2 ) 2 - c{£ x - £ 2 ) 2 ] , (3.6) 

where J\f is a normalization factor, 

M = [ W 2 (3 2 (a + b-W + c-W 2 )}~ 1 , (3.7) 

with W = \/s ' ( m N ~ m w)/^ m N^ P = (1 — ^m 2 N / s) 1 / 2 . The coefficients a, 6, c depend on 
the relative CP of the A^A^ system, and take the values 



+ 2 2 / 2 2 \2/o ( m AT + 2m^/) 2 

a + = m N m w (m N - m w ) (2 s — — ) , 

2m% 



b + = y/s m 3 N (m 2 N - 2m 2 v ) 2 (m 2 N - m 2 w ) , (r) C p = +1) 

c + = -m%(m 2 N - 2m 2 w f , (3.8) 

2 

aT = 2 m 2 w {m 2 N - m 2 v ) 2 (s(s - 2m 2 N ) (m 2 N + 2m 2 v ) 2 ) , 

b~ = y/s m N (m 2 N - 2m 2 v ) 2 (m 2 N - m 2 w )(s - 2m 2 N ) , (r) C p = -1) 

c~ = -m%(m 2 N - 2m 2 y ) 2 (s - 2m^) . (3.9) 

Notice that the ratios b + /a + and b~ /a~ are unequal (likewise the ratios c + /a + and c~/a~), 
although b + /c + = b~ /c~ . Thus the energy correlation of the two electrons in the final state is 



different for the cases r\cp = ±1. This is illustrated in Fig. 3 for the hypothetical parameters 
m N = 500 GeV, y/s = 1200 GeV. It may be noted that the factor f 2 = AcyCA does not 
appear in the spectrum (da / d£\d£ 2 ) , so that the energy correlation of e + e + dileptons is the 
same as that of e~e~. In the limit f3 — > 1 the term dominates and the '+' and '— ' cases 
are no more distinguishable. 



B. Angular Correlation 



Eq. (3.4) also allows a calculation of the correlated angular distribution of the final state 
e~e~ system. Defining the angles 6 1}2 as in Eq. (3.5), and integrating over all other variables, 
we find 



with 



T(m 2 N + 2m 2 w f s • /C? 1 /^ 2 



b F~i H~) ± ~ P ' / dcos °n { h ■ v - in \v) s ' 

a cos Bid cos 9 2 J I L 

+ (m 2 N - 2m 2 w f s ■ ( ±K e 2 1 K e 2 2 (cos 9 n (3 - cos 9 1 )(cos9 n (3 + cos0 2 ) 
+ /C^/C? 2 (l + /3cos0 n cos0i) + /C? 1 /C 2 ?2 (1 - 13 cos 9 n cos 9 2 ) ) 
2m 2 w s 2 • /C^/C? 2 (1 + cos 2 9 n p 2 ) - Am 2 N (m 2 N - 2m^) 2 



+ 



— cos #1 cos 9 2 ) 
+f 2 ■ 2(m% - Am 4 w ) s ■ 



COS( 



^(/CM 2 -/CM 2 ) : + 
K, e 2 1 )C e 1 2 cos9 1 + K 6 1 1 K e 2 2 cos9 2 : 



1(2) 



2Ai 



(2) 



(2) 



^ 2 (2)) 3/2 ' 



/c: 



'1(2) 



2^1(2) + ^1(2) 



(2) 



R2 \5/2 



(3.10) 



(3.11) 



^1(2) = 1 - (+)P cos 9 n cos 6>i(2) , #i(2) = sin 6> n sin 6>i(2) • 

The correlation (3.10) has been evaluated for m N = 500 GeV and y/s = 1200 GeV (using 
fi — 1 + 4sin 2 6V + 8 sin 4 / 2 = 1 - 4sin 2 6V) and is plotted in Fig. 4. There is a 
clear difference between the cases r]cp = ±1. The angular correlation in (3.10) becomes 
particularly transparent near the threshold f3 — > 0, where we obtain the analytic results 



9 



J_ , da + _ 1 1 ^ - 2m 2 

a+ ' [ d cos M cos 2 ^° ~ 4 ' [ 2 ' /i + 2m w 



7 • [ 1 + M • ft • (cos B x + cos 2 ) ] , (3.12) 



d(J x- 1 r i , ( m N - 2m w) 2 a a 



a- K d cos 9 t d cos 9 2 ,f3 ^° ~ 4 1 K + 2m 2 H/ ) 2 

. /2 - 2m^) 

+ A- (^ + 2^) - (cosgl + cosg2)] - (3 - 13) 

Notice that the distribution in the variables cos#i and cos# 2 becomes flat in the case rjcp = 
+ 1 when /2//1 is neglected. By contrast, there remains a nontrivial correlation for t]cp — — 1, 
even in the absence of / 2 . As before, the above results for e~e~ hold for e + e + if one replaces 

IV. UNLIKE-SIGN DILEPTONS: THE REACTION e + e~ -► JV1JV2 W+VF e+e 

Proceeding as in Sec. 3, the matrix element for the reaction e + e~ — > A^A^ — > 
Py+W^-e+e- (Fig. 2) is 

1 1 



M m = iA j; Af - 



(ft — mf + imiTi q 2 — m 2 + ^2^2 
[ A 2 ^(fci) 7 P ^1 7^ ^2 7<x g(! - 75) v t2 (k 2 ) (4.1) 

+ Ai m 1 m 2 u tl (h) 7 P 7^ 7 CT ^(1 - 75) v t2 (k 2 ) ] e* x p 3 (k 3 )e* x a A (k 4 ) ■ 

The same final state, produced via a Dirac particle-antiparticle pair (e + e~ — > NN 
W + W~e + e~), has the amplitude 

1 1 



M d = iA jl Af - 



q[ — m{ + irriiTi q 2 — m 2 + im 2 r 2 



M&i) 7p^i 7*102 la ^(1 -7s) uta (92) c^OaK^) • (4.2) 

Summing (averaging) over final (initial) polarizations, and using the narrow-width approxi- 
mation for the JVi, iV 2 propagators, we obain the squared matrix elements given below: 

= ^7-^^V%i-<)^V%2-<)^ 



2 (s-m|) 2 m^r n m Jvy m N T n Wi wy m 
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fi ■ { T(m' 2 N - m 2 w f{m 2 N + 2m 2 w f ■ s ± A(m 2 N - 2m 2 w f 

s (hk 2 )(qiq 2 ) - s (k 1 q 2 )(k 2 qi) - (hk 2 ) (qiq) (q 2 q) + (kik 2 )(qil)(q2l) 
+ (hq 2 )(k 2 q)(qiq) - (hq 2 ) (k 2 l) (qj) + (hq) (k 2 qi) (q 2 q) - (hi) (k 2 q ± ) (q 2 l) 
- (hq)(k 2 q)(qiq 2 ) + (k 1 l)(k 2 l)(qiq 2 ) ± m 2 N ((hq)(k 2 q) - (k 1 l)(k 2 l)) 
- 2 (m 2 N - m^)(m^ - 2m 2 w ) 2 (hq)(q 2 q) - (hl)(q 2 l) + (k 2 q)(qiq) - (k 2 l)(qil) 
+ 2 (m 2 N + Am 2 w /m 2 N )(m 2 N - m 2 w f [ (qiq)(q 2 q) - (qil)(q 2 l) } 
-2f 2 ■ ^ (m 2 N - m 2 w )(m A N - 4m^) ( ±(k 1 q)(q 1 l) - (hq)(q 2 l) 
T(hl)(qiq) + (hl)(q 2 q) T (k 2 q)(q 2 l) + (k 2 q)(q 1 l) 
±(k 2 l)(q 2 q) - (fc 2 Z)(gig) ) 

+2(m 2 N + Am A w /m 2 N )(m 2 N - m 2 w ) 2 \{q iq ){q 2 l) - (q 2 q)(qil)] ) [ . (4.3) 



\A\' 



m 



2 \2 



m N T n 



5{ql-m 2 N ) 



m N T n 



5{ql-m 2 N ) 



m 



N 



m 



w 



(m 2 N - 2m 2 w ) 2 ■ (k iq )(k 2 q) - (hl^l) 



m 



+2^-(m 2 N - m 2 w ){m 2 N - 2m 2 w ) ■ {kiq){q 2 q) - (hl)(q 2 l) + (k 2 q)(qiq) - (k 2 l){qil) 



m 



w 



+A(m N - m w ) ■ (qiq)(q 2 q) - (qil)(q 2 l) 



+ f: 



2 • 



m 



N 



777 



W 



(m 2 N - 2m 2 w ) 2 ■ (q iq )(k 2 l) - (k 2 q)(k 1 l) 



m 



+2^ r (m N - m w )(m N - 2m w ) ■ (k 1 q)(q 2 l) - (k 2 q)(qil) + (qiq)(k 2 l) - (q 2 q)(kil) 



+4(m : 



m 



,2 \2 

w) 



(QiQ)(Q21) ~ (q2q)(qil) 



(4.4) 



In complete analogy with the discussion of like-sign leptons (Sec. 3), we derive from the 
above equations the correlation in the energies and angles of the final e + e~ state. 



A. Energy Correlation 

The distribution in the scaled energies €i, £2 has the quadratic form given in Eq. (3.6), 
where the coefficients in the Majorana cases '+' and '— ' and the Dirac case 'd' now have 
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the values 



a+ = -m 2 N (m 2 N - m 2 w ) 2 (-^(m% + 4m^) - 2(m 2 N + 2m 2 w ) 2 \ , 
b + = -y/s m%(m 2 N - 2m 2 w ) 2 (m 2 N - m 2 w ) , 

c+ = m%(m 2 N - 2m 2 v ) 2 , 

1 / 4 

a' = -m 2 N (m 2 N - m 2 w ) 2 (s(s - 2m 2 N )(l + 4^f) - 4(m 2 N + 2m 2 w )- 

b~ = -y/sm N (m 2 N - 2m 2 w ) 2 {m 2 N - m 2 w )(s - 2m 2 N ) , 

c~ = m%{m 2 N - 2m 2 v ) 2 (s - 2m 2 N ) , 

/ 4 



a = 



(m 2 N - m 2 w ) 2 ^4(s -m 2 N )(s- Am 2 N ) 



rn 



N 



~(m 2 N 



- 2m 2 v )(2(s - 4m^)m^ - m 2 N (m 2 N - 2m^))^ , 
b d = ^/s m N {m 2 N - 2m 2 v ){m 2 N - m^) (astu^ - (m^ + 14m^)m^ 
c rf = m 4 w K-2 m ^) 2 (s-3m^ . (4.5) 

The corresponding three distributions are plotted in Fig. 5. As in the case of like-sign 
dileptons, the e + e~ pairs have distinct correlations for rjcp = ±1. A comparison of the 
Majorana cases with the Dirac case reveals an interesting difference. In the Majorana cases 
the total e + e~ energy Y = E\ + £2 is distributed symmetrically around the mid-point of this 
variable Y = l/2(Y min + Y max ). By contrast, e + e~ pairs resulting from Dirac NN primary 
state have a total energy distribution that is unsymmetric around the mid-point. 

B. Angle Correlation 

In analogy to the distribution da/d cos 9\d cos 6 2 obtained for e~e~ pairs (Eq. 3.10), the 
result for unlike-sign dileptons e + e~ is 

( , , ) ± ~ P- I dcoseJh- 

d cos did cos V 2 J l 

' T2m 2 N {m 2 N + 2m 2 w ) 2 s ■ K%K% 

- 2m 2 N (m 2 N - 2m 2 w ) 2 s ■ ( ±K% K% (cos 6 n p - cos #1) (cos6 n (3 + cos9 2 ) 
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+ /C^/Cf (1 +/3cos0 n cos0i) +/C? 1 /C^ 2 (1 - f3 cos 9 n cos 9 2 ) ) 

+ (m% + 4m 4 w ) s 2 ■ K^ICl 2 (1 + cos 2 6 n {3 2 ) + 8m%(m 2 N - 2m 

K%K%(\ - cos 0i cos0 2 ) 

+/ 2 • 2(rr4 - 4m*,) s • [ -/CM 2 • s cos6 n (3 

cos 9 n f3 (/C^/C? 2 + /C? 1 /C 2 ?2 ) : + 
/C^ 1 /C? 2 cos 0i - /C? 1 /C^ 2 cos 2 : - 



2 \2 
14^7 



(4.6) 



h ^ ^ ~ /?•/ rfcos0 n {/i- 

dcos0idcos0 2 J I 

' s 2 m i w -K e l 1 K B l 2 ■ (l + cos0 n/ 3 2 ) 

2 -2m 2 w )s- (/C^ 1 /C^ 2 (l + /3cos0 n cos0i) + /C? 1 /C^ 2 (l-/?cos0 n cos0 2 ) 



, 2 2/2 

+ m N m w (m N 
+ m%(m 2 N - 2171^) -/v 2 

2 m^-/CM 2 -cos0 n /3 



^(m^ - 2m^) 2 • /C^ 2 • (1 - cos0i cos0 2 ) 



-f-2 ■ \ 2s 1 m% / - , | - ; „ . 

+ m^m^m^ - 2777^) s • ( /C^/C^ 2 (cos0 n/ 9 + cos0i) + /C? 1 /C 2 ,2 (cos0„/3 - cos0 2 ) ) 

+ m^(m 2 v -2m 2 y ) 2 -/C^ 1 /C^ 2 -(cos0i-cos0 2 )] J . (4.7) 

As usual, the indices '+', '— ' and 'd' differentiate between the Majorana cases r/cp = +1, 
— 1 and the Dirac case. The angle-correlations expressed by Eqs. (4.6), (4.7) are plotted 
in Fig. 6., where the differences between the three cases are obvious. Close to threshold 
{(5 — > 0), the correlation between cos0i and cos0 2 can be presented in analytic form 

da 11*. ™2 o_2 



a+ V rfcos0irfcos0 2 ;/3 



\+ 1 r t . 1 /2 — 2777^, 

l «"r |1+ 2i'i5rt4' M, " 0S<,l)l ' 



(4.8) 



da 



-) d - — 



rT rf ' ^ d COS 0id COS 2 (7- " V d COS 01 d COS 2 



(4.9) 



1 (m^ - 2m 2 y ) 2 / 2 (m^ - 2771^) 

4 ( m % + 2<)» ' C0S * C0S * + A ' KTK) ' (C0S <' 1 " ^ 1 ' 

In this limit, the cases '+' and '— ' remain distinct, but the case r\ C p = —1 converges to the 
Dirac case. 
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V. COMMENTS 

We comment briefly on some other papers which have a partial overlap with the consid- 
erations presented above. 

(i) Our discussion of the production reaction e + e~ — > NiN 2 follows very closely that given 

in Ref. Our results for da/dfl given in the Appendix (Eqs. (A.1)-(A.6)) essentially 
coincide with those in this paper, with two minor differences: The angular distributions 
for the case of two distinct Majorana particles with the same CP-parity, as well as 
for the case of two distinct Dirac particles (Eq. (4E) and (4D) in Ref. Q), are slightly 
different from our distributions, presented in the Appendix (Eq. (A. 2) and (A. 3)). 

(ii) The cross section for the Majorana process e + e~ — *> NiN 2 , with mi = m 2 and rjcp = +1 

calculated in Ref. [[| agrees with that obtained in this paper. However the Dirac case 
e + e~ — > NN (Eq. (2) of Ref. ||) differs from our result (Eq. (A.5)), as also noted in 
Ref. 1. 

(iii) The spin-summed differential cross section for the Majorana process e + e~ — > NiN 2 
(with mi = m 2 , r/cp = +1) calculated in the present paper, as well as in Refs. 
differs from that given in Ref. 0, but agrees with the results given in Refs. P,|5|,|ID[. 

(iv) Our analysis of heavy Majorana production and decay has been essentially model- 
independent. Discussions in the context of specific gauge models, based on SU{2) L x 
SU{2)r x U{1) or E(6) symmetries, may be found in Refs. [P. |ITi| . |TT|1 . 



14 



APPENDIX A: DIFFERENTIAL CROSS SECTION FOR e+e~ -> N ± N 2 

Following Ref. Q, we consider the following five cases, in which N% and N 2 are 
A distinct Dirac particles. 

B distinct Majorana particles with the same CP-parity. 
C distinct Majorana particles with opposite CP-parity. 
D Dirac particle-antiparticle pair. 
E identical Majorana particles. 

Choosing the Ni direction in the e + e~ cm. system to be the z-axis, and the e~-beam 
direction to be at an angle 6 (= scattering angle), the momenta (gi, q 2 ) and spins (ti, t 2 ) of 
Ni and N 2 have components 

N i : qt = (t^IjOj 0,7/3mi) , 

= (j(3n z ,n x ,n y ,jn z ) , 
N 2 : ^ = ( 7 'm 2 ,0,0,-7 , /3 / m 2 ) , 

The differential cross sections are (with (3 — (1 — 4m 2 /,s) 1 / 2 , /?' = (1 — 4m 2 /,s) 1//2 , 7 = 
(1 - /3 2 )- 1 / 2 , 7' = (1 - {3' 2 )- 1 ' 2 , X(x, y,z) = [x 2 + y 2 + z 2 - 2(xy + yz + zx) J 1 / 2 ) 

= f^f l^)| 2 [1 - (ml - ml) 2 Is 2 ] \(s,mlm 2 2 ) 

{A [ (1 + C 2 (3'(3) - ((3n' z + (3'n z )(l + C 2 ) - ((3'nJ^ + (3n'Ji)SC 

+ n z n' z (C 2 + (3(3') + (n x n'J 1 + n' x n z /i)SC + n x n' x S 2 } 
+ h [ C((3 + (3') - (n z + <)C(1 + (30) - [n x j 1 + n' x /j')S 
+ n z n' z C((3 + (3') + S{P'n x n' z h + [3n' x n z /j') } } , (A.2) 
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{%)b = ^| \R(s)\ 2 [1 - (ml - m 2 2 ) 2 /s 2 } X(s,ml,ml) 



{ /i [ nX^ 2 (l/77' " 1) + + n,n' 2 (^' - C 2 (l/ 7 7' - 1)) 

+ (n x n' z - n' x n z )SC(lh ~ 1/V) + C 2 (3(3' - 1/77' + 1 ] 
+ h [ K " <)S(l/i - 1/7) + (n, + <)C(l/ 77 ' - (3(3' - 1) ] } , (A.3) 



= w |i?(s)|2 [1 " {ml ~ mlf,s2] X ^< m2 ^ 

fi [ n x n' x S 2 (l/^' + 1) - n y n' y S 2 {3{3> + n z n' z ((3(3' + C 2 (l/ 7 7' + 1)) 

+ (n x n' z + n' x n x )SC(lh + W) + + 1/77' + 1 ] 

h [ (n x + <)5(l/ 7 + 1/7') + (n, + <)C(l/ 7 7' + (3(3' + 1) ] } , (A.4) 



{ /1 [(1 + C 2 /3 2 ) - (n 2 + + C 2 ) - (n x + n' x )SCPh + n z n' z {C 2 + /3 2 ) 

+ (nX + n 2 <)5C/ 7 + nX5 2 /7 2 ] + h [ 2C (3 - (n z + n' z )C (1 + (3 2 ) 

- (n x + n' x )Sft + 2n z n' z C(3 + (n x n' z + n' x n z )S(3h ] J , (A.5) 



/1 [{n y n'-n x n' x )S 2 + {l + n z n' z ){l + C 2 )} - f 2 2(n z + n' z )C \ . (A.6) 



APPENDIX B: MATRIX ELEMENTS FOR e+e -► A^iVs -► 6^^^+ 

1. Like Sign Dileptons 

The matrix element for the reaction (Fig. 2) 

e + (p 2 ,t 2 ) + e-(p 1 ,t 1 ) -> e'(A; 1 ) e -(A; 2 )^ + (A;3,A3)^ + (A;4,A4) (B.l) 



is 
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+ ^2 1 

,I-mI + im 2 r 2 ^2 (1 + 75) ^ (A;2) 
- Ai u t2 (k 2 ) la hi - In) o ^ 2 + , m L f 7* J(l - 7s) 



2 



-^i + m 1 1 



qf — ml + imiTi 
Upon rearrangement, this gives the matrix element in Eq. (3.2) 



1 P 2(1 + 75) t*(*i) } e£(*3)<£(*4) . (B.2) 



2. Unlike Sign Dileptons 



The matrix element for the Majorana- mediated process (Fig. 2) 



e + ( P 2,t2) + e~(p 1 ,t l ) -> e-(A; 1 ) e -(A; 2 )VU + (A;3,A3)lU + (A; 4 ,A4) (B.3) 



is 



* m = -a ^Ar{A 2 , tl (M 7^(1 -75) ^^7^(1-7,) 

-^ 2 + m 2 1 
-2 — 9 1 -p 7^(1- 75) v t2 (k 2 ) 



g 2 — m 2 + im 2 T 2 2 

" Al ^ {h) ^ V 1 - 7s) g -t^Ur, 7 ' ^ (1 + 75) 

2 "latTlLr. 7. £(1 - 75) t/Jfca) j cS(*3)^(*d) • (B.4) 



g 2 - mi + im 2 r 2 2 
Upon rearrangement, this gives the matrix element in Eq. (4.1). 
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FIGURE CAPTIONS 

1. Feynman diagram for the reaction e + e~ NiN 2 

2. Diagram showing the sequential process e + e~ — > NiN 2 — > e ± e~W Zf W + 

3. Energy correlation of the e~e~ lepton pair in the reaction e + e~ — > NiN 2 — > 
e~W + e~W + , for the cases (a) r^p = +1, (b) ?y C p = —1. (Paramers for this and 
succeeding figures: y/s = 1.2 TeV, rrijy = 500 GeV.) 

4. Angle correlation of e~e~ dileptons in e + e~ — > A^A^ — * e~W + e~W + , for (a) ??cp = 
+ 1, (b) ?fcp = -1 

5. Energy correlation of e~e + dileptons in e + e~ N\N 2 — > e~W + e + W~: (a) Majorana 
pair, r/cp = +1, (b) Majorana pair, rjcp = — 1, (c) Dirac A/W-pair 

6. Angle correlation of e~e + dileptons in e + e~ — > A1A2 — > e~W + e + W~: (a) Majorana 
pair, 7/cp = +1, (b) Majorana pair, r/cp = —1, (c) Dirac A^A'-pair 
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FIGURES 
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FIG. 3. 
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e e Final State: Angular Correlation 
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FIG. 4. 
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e + e Final State: Energy Correlation 
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FIG. 5. 
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e + e Final State: Angular Correlation 



Majorana (+) Majorana (— ) 




Dirac 




FIG. 6. 
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